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a^ : ABSTRACT 

o 

^ I We present a method for extracting the mass of the top quark from the di-lepton 

decays of top anti-top quark pairs. In this decay channel two neutrinos remain un- 
detected. Extraction of the mass of the top quark by kinematic reconstruction is not 
possible because the event is under-constrained. We therefore employ a dynamical 



D I likelihood method to solve the problem. 

> 

^ ' 1 The di-Iepton event topology 

In the di-lepton channel top anti-top quark pairs decay via t Wb, followed by 
W lu. We can identify and measure the 4-momenta of the jets^ and the charged 
leptons. The two neutrinos in the event remain undetected. The vector sum of their 
transverse momenta can be inferred from the observed missing pt- This leaves us 
with 14 observables {a} out of 18 values {v} needed to describe the six particle final 
state. In order to constrain the ti kinematics, three additional constraints on the 
invariant masses of final state particle combinations are introduced: 
-k m{l^,vi^) = m{l2,i'i^) = mw, 



^Both jets are assumed to be 6-jets 
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In the end the set of equations are still under-constrained by one equation. Ex- 
traction of the mass of the top quark by kinematic reconstruction is therefore not 
possible. 

2 The Analysis 

This analysis was developed during Run I[T]. Ideally we would like to calculate the 
probabihty to measure the 14 observables {o}, given the top quark mass mt- 

P{{o}\mt)oc f f{x)f{x)\M\^p{{o}\{v})d^^{v}dxdx (1) 

J{v} 

is the probability density function {pdf) for measuring {o} given the 18 final state 
parameters {v} for a hypothetical rrit and for the parton (anti-parton) momentum 
fraction x (x). Here f{x) {f{x)) is the proton (anti-proton) parton distribution 
function and A4 is the matrix element for the process. A hypothetical value of nit 
is used as the last constraint. We then obtain up to 4 real solutions for the u and 
u momentajS]. There is a two- fold combinatoric ambiguity in pairing the charged 
leptons and 6-jets. As a result up to eight solutions of the neutrino momenta are 
possible. Instead of rigorously computing the pdf in Equation [H event weightsj^jH] 
are used to characterize the physics. In this poster we calibrate the effect of this 
simplification by comparing the measured mass of the top quark obtained from the 
likelihood fits to simulated events, versus the value used in their generation. 

For every event, a weight (Wfc) which is a function of the hypothesized 
mass of the top quark {rrit), is derived [2] |S] that corresponds to the k^'^ neutrino 
momentum solution. The event weight 

Wk{mt) oc f{x)f{x)p{E'\mt)p{E'\mt), (2) 

represents the likelihood to observe the event for a given mt- Here p{E'\mt) is the 
probability density function for the energy of the charged lepton to be E' in the rest 
frame of the top quark of mass m^P]. The quantity p{E'\'mt) is the analogue for the 
anti-top quark. With n solutions for the neutrino momenta the total event weight 
W is defined as: 

n 

W = normalization ^ Wfe (3) 

k=l 

The value of rrit at which the weight curve peaks is used as the mass 
estimator in the analysis. All simulated events are first filtered using the same 
criteria as for data events|2|. Figure ^left) shows the weight curve for a simulated 
event with input mt = 175.0 GeV. The peak values from many events generated 
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Figure 1: (left) Example weight distribution from a simulated event generated with 
rrit = 175.0 GeV. (right) Template constructed from peak values obtained from 
many events each having input rut = 175.0 GeV. 

with the same rrit are binned into histograms to construct templates, one of which 
is shown in Figure fright). 

It is not possible to generate MC events with continuously varying input 
rrit. Seven distinct rrit values are used to generate the signal MC template distri- 
butions from 120.0 GeV to 230.0 GeV. Templates representing contamination from 
background processes are also constructed and then added to the signal templates 
in the proportion of the expected background Normalized templates represent 
the likelihood of observing the particular event if the mass of the top quark equals 
nit. 

Sets of eight MC events are used to construct simulated experiments^. A 
binned maximum likelihood fit is performed for each such ensemble using the tem- 
plate distributions. One such fit to ensemble events is illustrated in Figure Efleft). 
Many different ensembles are constructed, and the exercise is repeated many times to 
obtain a distribution of maximum likelihood estimates (MLEs) as in FigureEf center). 
The mean of this distribution yields the mean fitted estimate of rrit. Repeating the 
same exercise for different input MC rrit yields values which can be used to check 
the performance of the method. A straight line parameterized as: 



gives the best fit to the ensemble test results for = 1.00±0.04 andpo = 175.8±0.08 
GeV. This fit to the set of points, shown in Figure |2fright), is consistent with a 
straight line of unit slope and 175.0 GeV offset. 

^Eight emu events observed from 142.73 pb~^ data|S]. 



fitted mass = pi ■ (input mass — 175.0 GeV) + po GeV 
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Figure 2: (left) Binned maximum log-likelihood fit of an ensemble of events to the 
template histograms yield the most likely estimate (MLE) of rrit for that ensemble, 
(center) A histogram of MLEs derived from 100 different ensembles, all signal events 
used have an input rrit = 175.0 GeV. (right) Calibration of a series of different input 
nit values. The 175 GeV input MC value follows from the histogram in the center. 

3 Conclusion 

The calibration of this method verifies its performance, and indicates absence of bias 
in the method. 
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